Many-body interband tunneling as a witness for complex dynamics in the 

Bose-Hubbard model 
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A perturbative model is studied for the tunneling of many-particle states from the ground band 
to the first excited energy band, mimicking Landau-Zener decay for ultracold, spinless atoms in 
quasi-one dimensional optical lattices subjected to a tunable tilting force. The distributions of the 
computed tunneling rates provide an independent and experimentally accessible signature of the 
regular-chaotic transition in the strongly correlated many-body dynamics of the ground band. 
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The experimental advances in atom and quantum op- 
tics allow the experimentalist to directly study a plethora 
of minimal models which have been developed to de- 
scribe usually much more complex phenomena occurring 
in solid states P, 0, S| • Bose-Einstein condensates loaded 
into optical lattices, which perfectly realize spatially peri- 
odic potentials, are used, e.g., to implement the Wannier- 
Stark problem ass as a paradigm of quantum trans- 
port where atoms move in a tilted lattice. Up to now all 
experiments on the Wannier-Stark system with ultracold 
atoms have been performed in a regime where atom-atom 
interactions are either negligible [J] or reduce to an effec- 
tive mean- field description [5j, 0] • State-of-the-art setups 
are, however, capable to achieve small filling factors of 
the order of one atom per lattice site Q. Moreover, the 
atom-atom interactions can be tuned by the transversal 
confinement and by Feshbach resonances [H, 0] , resulting 
in strong interaction-induced correlations. 

The regime of strong correlations in the Wannier-Stark 
system was addressed in 0, revealing the sensitive 
dependence of the system's dynamics on the Stark force 
F. The single-band Bose-Hubbard model of [tj [l(| is 
defined by the following Hamiltonian with the creation 
a] annihilation a^i, and number operators hi t i for the 
first band of a lattice I = 1 . . . L: 
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FIG. 1: (a,b) CDF (stairs) and P(s) (stairs in insets) for N = 
5 atoms, L — 8, lattice depth V = 10 recoil energies (fixing 
Ji = 0.038), Ui = 0.032, F ~ 0.063 (a) and 0.021 (b), with 
WD (solid) and Poisson distributions (dashed), (c) \ 2 test 
with values close to zero for good WD statistics. The dashed 
line marks the transition to quantum chaos as F is tuned, (d) 
variance of the number of levels in intervals of length dE (with 
normalized mean spacing), for the cases of (a) (squares) and 
(b) (circles), with the random matrix predictions for Poisson 
(dashed) and WD (solid) [H. 



X 



— ni t i (n /4 - 1) 



(1) 

A transition from a regular dynamical (dominated by F) 
to a quantum chaotic regime (with comparable values of 
Jx,U\ , F) was found 0, [n| ■ The transition was quan- 
titatively studied using the distribution of the spacings 
between next nearest eigenenergies of the Hamiltonian 
(fl]) . This analysis [13] verifies that the normalized level 
spacings s = AE/AE obey a Poisson (P(s) = exp(— s) ) 
and a Wigner-Dyson (WD: P(s) = s-k/2 exp(— its 2 / 4) ) 
distribution in the regular and chaotic case, respectively 
[ill ]. P( s ) an d the cumulative distribution functions 
(CDF: C(s) ee J o s ds' P(s')) are shown for typical cases m 
Fig. [TJ where we scanned F to emphasize the crossover 



between the regular and the chaotic regime. Statistical 
tests are also shown which confirm the analysis of [^, [13] 
in a more systematical manner 12 1 . 

As shown in [§] , the strong correlations in the quantum 
chaotic regime induce a fast and irreversible decay of the 
Bloch oscillations, which otherwise would persist in the 
ideal, non-interacting case. Therefore, the crossover be- 
tween the two regimes discussed above could be measured 
in experiments by observing just the mean momentum as 
a function of time. Here we introduce a new, robust and 
hence also experimentally accessible prediction for this 
crossover. In the presence of strong interactions parame- 
terized by Ui , the single-band model should be extended 
to allow for interband transitions [13] , as recently realized 
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at F = in experiments with fermionic interacting atoms 
Q . Instead of using a numerically hardly tractable com- 
plete many-bands model, we introduce a perturbative de- 
cay of the many-particles modes in the ground band to 
a second energy band. Our novel approach to study the 
Landau-Zener-like tunne ling between the first and the 
second band [H, d, 0, Ei , 15 1 leads to predictions for the 
expected decay rates and their statistical distributions. 
As we will show, the latter are drastically affected by the 
dynamics in the ground band, and they therefore provide 
a measurable witness for the regular-chaotic transition. 

We first derive the individual decay rates of the dom- 
inating interband coupling channels. These decay rates 
will serve to effectively open the single-band model ([1]) 
for mimicking losses arising from the interband coupling. 
Our analysis starts from the following "unperturbed" 
Hamiltonian for the first two bands: 
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£i his + £2 nj,2 - 4f-(ai+i,2 f «;,2 + h.c.) 
+Fl{h h i + h lt2 ) + \h l<l (h l!l -l)] . (2) 



For a moment, we neglect the hopping in the lower band, 
where the single-particle Wannier functions [3] are more 
localized than in the upper band. In the latter we neglect 
the interactions, since initially only a few particles pop- 
ulate the excited levels. A closer analysis of the full two- 
bands system [l2| shows that there are two dominating 
mechanisms that promote particles to the second band. 
The first one is a single-particle dipole coupling arising 
from the force term: 



Hi = F ■ D ^2 



a;,2 f a,i + a/V d L2 



Here, \w) = Y^tZ-oa^rn-wiVh]/ F)d m . 2 ^\va,c) is the 
single-particle eigenstate for the Wannier-Stark problem, 
localized around the site w in the second band, with the 
Bessel function of the first kind J m {x) [3|. 

The expectation value of ((3]) for \b; N) of the first band, 
equal to the first-order SE(b), is zero because the opera- 
tor does not conserve the number of particles within the 
bands. The decay width at first-order is given by the ma- 
trix element of the perturbation between the initial and 
final state according to Fermi's Golden Rule, and only 
the first term in ([3]) gives a nonzero contribution [H3 |: 
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The <$(•,•) functions act as a selection rule for the Fock 
states that are coupled by the perturbation. The tun- 
neling mechanism does not include any income of en- 
ergy from an external source, so the initial and final 
energies E (b) = (vac| (b\H Q \b) |vac) and E (b',w) = 
(w\(b'\H \b'}\w}, respectively, must be equal as required 
by the Golden Rule. The condition on the energy conser- 
vation is, however, relaxed to account for the uncertainty 
AE(b) of the unperturbed energy levels of the initial and 
final states in the lower band arising from the hopping in 
this band initially neglected in ^j. A detailed derivation 
is given in (l2| . and here we only state the result: 

AE(b) = 2vr ( Ji /2) 2 ^ AE(b -> b') = 

b< 

2tt(Ji/2) 2 rif 5(ni+Ai + l,ni). (7) 

l Al=±l 



where D depends only on the lattice depth V (measured 
in recoil energies according to the definition in Q). The 
second one is a many-body effect, describing two particles 
of the first band entering the second band together: 



H 2 = ^ {^ai^ ai^ a l} i a lA + (!<-> 2)) 



(4) 
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The cross-band interaction is characterized by the pa- 
rameter U x = a s J dxxlxl - °- 5C/ i ( for V = 3...10) 
for Ui= a s J dx\ii with renormalized scattering 
length d s [3, U3] and the Wannier functions %i,2 local- 
ized in each well for the first or second band. To justify 
the following perturbative approach, it is crucial to real- 
ize that the terms ((31) and Q must be small compared 
with the band gap A = e 2 — £i (not necessarily small 
with respect to the single band terms in (TTJ), and indeed 
FD, U x ,Ui <C A for the parameters considered here. 

For the first perturbation, the decay channel of a given 
unperturbed Fock state labelled \b) (with a total number 
of atoms N and nn atoms in an arbitrary well h) is 

|6; N) $ |vac) -> \b'; N-l)®\w), n' h =n h -l. (5) 



The level density p(E, b) around the unperturbed en- 
ergy Eo(b) of a Fock state \b) is then approximated by 
a rectangular profile, of width AE(b) and unit area: 
p(E,b) = x{\E-E (b)\ < AE(b)/2} / AE{b). The re- 
laxed energy conservation rule selects from §5§ the set K 
of permitted decay channels (h, w) parameterized by the 
two indices h.w such that: 



E (b',w)-E (b) 



F(h -w)-Ui (n h - 1) 



AE(b)+AE{b' 
2 



AE{b)+AE(b') 
2 



(8) 



Hence the energy A required to promote a particle to the 
second band is supplied by the decrease of the interaction 
(cx Ui ) and by the work of the force (ex F) exerted on 
the promoted particle. 

The total width Ti(b) for the decay via the allowed 
channels K, is proportional to the square of the matrix 
element and to the level density p(E,b): 



Ti{b) = 2^{FDfY. 



(h,w)eK 



AE(b)AE( 



AE(b') } 



(9) 
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J m {x) significantly contributes only for |m| 
Ui , AE(b) <C A, the energy conservation is roughly 
given by |A| ~ F(h — w). Requiring that the Bessel 
function in §§§ is substantially larger than zero, we ob- 
tain the inequality |A| < | «7a ] - The last condition does 
not depend on F, since a twofold effect is at work: a 
stronger force produces a larger energy gain when a par- 
ticle moves along the lattice, but the extension | J2/-FI of 
the single-particle state shrinks. Therefore, increasing F 
results in an increased energy matching and a strongly 
reduced "geometrical" matching. For 3 < V < 26, we 
have |A| — I J 2 | > 1.0 such that the energy matching 
cannot be realized by just tuning the lattice depth. The 
decay can, however, be activated by an increase of the in- 
teractions, which can be experimentally achieved by act- 
ing on the transversal confining potential of a quasi-one 
dimensional lattice, or by a Feshbach resonance [8|. In 
the calculations presented below, we augmented Ui used 
in |T3| by a factor of order 10, and as noted in the in- 
troduction, a similar increase of the interaction strength 
was used in the experiment to promote fcrmions to higher 
bands 0] , in close analogy to the here described field- and 
interaction-induced interband coupling of bosons. 

The second term (j4} is treated in a similar way, with 
the difference that two particles are promoted to the sec- 
ond band, and the position of the second single-particle 
state \w') is an additional degree of freedom for the tran- 
sition. The decay channels are: 

|6, N) ® |vac) -f \b', N - 2) <g> \w, w') ; n' h =n h -2. (10) 

The energy matching selects a set K of decay channels, 
parameterized by the three site indices h, w, w': 
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(h,w,w') S K s.t. E (b',w,w')-E (b) = 
= 2A - F(2h -w-w')-Ui (2n h - 3) 
AE(b) + AE(b') AE(b) + AE{b') 



(11) 



The computation of the matrix element yields [1^ 
2 



r 2 (6) = 2tt 



E 



(h,w,w')£K 



^ yJ h — w ( p ' ) 



J k-w' ( p ) • 4n/i (n h 1) A1 . , jA Li , , 



(12) 



With respect to ([9]), the additional degree of freedom 
w' results in a summation over all possible values of 
w — w' . This follows from the possibility to conserve the 
energy even if a particle is pushed far, if the other parti- 
cle is pushed almost equally far in the opposite direction. 
Since the decay widths in (fT2|) depend on the product of 
two (rapidly decaying) Bessel functions - again a "geo- 
metrical" matching condition - we apply the truncation 
\w — w'\ < I J2/FI, to reduce the formula to a finite form. 

We can now compute the total width Tp(b) = Ti(&) + 
^2(0) defined by the two analyzed coupling processes for 
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FIG. 2: (a,c,e) CDF from Re{Sj} (stairs), together with 
WD (solid) and Poisson predictions (dashed). (b,d,f) distri- 
butions of the logarithm of the rates. In (a,b), (c,d), (e,f) 
F ~ 0.17,0.31,0.47, respectively, with (N, L) = (7,6),F = 
3, Ui — 0.2 (fixing U x — 0.1). In the regular regime (f), 
a log-normal distribution (dotted) well fits the data, with a 
scaling P(F) oc T~ x for the largest F (dashed line in the inset 
of (f) with x = 1). In the chaotic case, a global power-law 
behavior with x ~ 2 is found (dashed line in the inset of (b)). 



each basis state |6) of the single-band problem given in 
(PJ. The Tp(b) are inserted as complex potentials in the 
diagonal of the single-band Hamiltonian matrix. After 
a gauge transform that recovers the translational invari- 
ance of the problem (see [l(], E2| for details), the latter 
matrix is used to compute the evolution operator over one 
Bloch period Tb, which is finally diagonalized to obtain 
its eigenphases exp (— iEjT-o)- Along with the statistics 
of the level spacings defined by Re {Ej}, the Figs. H]and[3] 
analyze the statistical distributions of the tunneling rates 
Tj = — 21m for some paradigmatic cases. All rates 
are much smaller than unity, which a posteriori is fully 
consistent with our perturbative approach. 

To observe what happens at the regular-chaotic tran- 
sition (c.f. Fig. Q]), we scan F in Fig. [21 and as F in- 
creases, the average decay increases by orders of mag- 
nitude, while the distributions broaden. The large in- 
crease of the rates is due to an improved energy match- 
ing, when F supplies the necessary energy to promote 
particles to the second band. For the parameters of 
Fig-Hi the single-particle Landau-Zener formula 14| gives 
r LZ - F/(27r)exp[-7r 2 A 2 /(8F)] - lO" 23 , 10" 12 , 10~ 8 
for (b,d,f). This huge variation, typical of semiclassical 
formulae, implies that there are possibly parameters for 
which our results are comparable to the single-particle 
prediction, but, in general, the many-particle effects can- 
not be neglected. Moreover, mean-field treatments of 
the Landau-Zener tunneling at best predict a shift of T 
0)0] i but cannot account for their distributions. 

In the chaotic regime, the Fock states are strongly 
mixed by the dynamics J^, [ljl [13] and a fast decaying 
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FIG. 3: (a,c) rate distributions in the chaotic regime with F ~ 
0.17, Ui = 0.2 (U x — 0.1), together with the corresponding 
unsealed P(T) in (b,d). In (a,b) (N,L) = (7,6),V = 4, and 
in (c,d) (N, L) = (9,8),V = 3. Power-laws P(r) oc T~ x are 
found with x ~ 2 (dashed lines in (b,d)). 



In summary, our perturbative opening of the single- 
band Wannier-Stark system allows one to study Landau- 
Zener-like interband tunneling within a many-body de- 
scription of the dynamics of ultracold atoms. The statis- 
tical characterization of the tunneling rates (mean values 
and form of the distributions) provides clear and robust 
signatures of the regular-to-chaotic transition for future 
experiments. A more detailed analysis of the interband 
coupling in a full-blown model, in which at least two 
bands are completely included, calls for huge computa- 
tional resources to access the complete quantum spec- 
tra. Nonetheless, our results are a first step in the di- 
rection of studies for which "horizontal" and "vertical" 
quantum transport along the lattice are simultaneously 
present and influence each other in a complex manner. 
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Fock state can act as a privileged decay channel for many 
eigenstates. Many states then share similar rates, lead- 
ing to thinner distributions. Therefore, the thinner dis- 
tribution of Fig. [2] (b) is a direct signature of the chaotic 
dynamics evidenced in (a), as compared with the regular 
case in (e,f). In Fig. [2] (f), we found a good agreement 
with the expected log-normal distribution of decay rates 
[ill (or of the similarly behaving conductance [l7|) in 
the regular regime. There the system shows nearly per- 
fect Bloch oscillations [9], and the motion of the atoms 
is localized along the lattice 
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We can even detect a 
qualitative crossover to a power-law P(r) oc T^ 1 in the 
right tail of the distribution, as predicted from localiza- 
tion theory [H, [H The distributions in Figs. [2] (b) 
and follow the expected power-law for open quantum 
chaotic systems in the diffusive regime [18J. The expo- 
nents x are, however, nonuniversal and depend on the 
opening of the system. In our case, the decay channels 
are defined by the interband coupling, which in a sense 
attaches "leads" to all lattice sites within the sample. 
Going along with the regular-to-chaotic transition in the 
lower band of our model (from Fig. [2] (f) to (b), or to 
Fig. [3]) the Y distributions transform from a log-normal 
to a power-law with x ~ 2, in close analogy to the tran- 
sition from Anderson-localized to diffusive dynamics in 
open disordered systems 18, 20j. 
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